Perfect mirror transport protocol with higher dimensional quantum chains 
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A globally controlled scheme for quantum transport is proposed. The scheme works on a ID chain of nearest 
neighbor coupled systems of qudits (finite dimension), or qunats (continuous variable), taking any arbitrary 
initial quantum state of the chain and producing a final quantum state which is perfectly spatially mirrored about 
the mid-point of the chain. As a particular novel application, the method can be used to transport continuous 
variable quantum states. A physical realization is proposed where it is shown how the quantum states of the 
microwave fields held in a chain of driven superconducting coplanar waveguides can experience quantum mirror 
transport when coupled by switchable Cooper Pair Boxes. 
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A fundamental issue in quantum information theory and the 
development of a quantum computers is the transfer of a quan- 
tum state through a reliable, high fidelity, method. Such meth- 
ods are known as quantum wires, or q-wires. Q-wire design 
should avoid the detailed control of its components as this 
may introduce additional errors which corrupt the transport 
and typically implies challenges in implementation. Ideally 
such quantum transport should be capable of transporting en- 
tire entangled registers of quantum information with near per- 
fect fidelity and operate with the minimum of external control. 
Recently interest in q-wires has grown considerably due to the 
advent of higher dimensional quantum systems (qudits or C V, 
continuous variables) as components in quantum processors. 

So far, the only proposals for a mirror transport involved 
qubit states [ 1 1 and generally required significant resources to 
achieve useful transport fidelity |j2j O. In this Letter a new 
protocol is described that for the first time formulates and 
solves the problem of quantum mirroring in chains consisting 
of nearest neighbor coupled identical qudit systems or infinite 
dimensional qunat [Continuous Variable (CV)] systems. It is 
based on a cellular automata-like scheme for perfect mirror 
transport J4], spatially inverting a quantum state distributed 
along a chain via the application of homogenous single qubit 
pulses and Ising-like nearest neighbor interactions. 

Qudits and qunats often lead to features which have cer- 
tain advantages over qubits |6). Fault tolerant quantum 
computation has been shown for qudits [7|, and schemes for 
trapped ion qudit quantum computation have been developed 
J8). One-way quantum computing via qudit cluster states 
is possible, as is the formulation of qudit decoherence free 
subspaces and topologically protected qudit memories (9)- 
Moreover, continuous variable (CV) quantum computation 
has been widely studied in recent years [ 10], while CV quan- 
tum transport has also attracted recent interest ifTTl . Our CV 
results below provide an alternative method for perfect quan- 
tum transport in CV chains. 

The aim is to develop a quantum mirror transport protocol 
for qudit and CV systems. The qubit case was considered 
previously [4 |. In this Letter a scheme suitable for qudits and 
CVs is developed. Physical realization is proposed involving 
perfect spatial mirroring of continuous variables in Circuit- 



QED 1 12], spatially mirroring a quantum state of microwave 
photons trapped in a chain of coplanar waveguides (CPWs), 
coupled via Cooper pair boxes (CPBs). 

Qudit mirror transport. It will prove useful to recall some 
basic properties of qudit quantum gates lfT3l [141 . A qudit is 
a state of a d— level system described by a d— dimensional 
Hilbert space fid with basis {\k); k = 0, 1, • • • , d— 1}. Gen- 
eralized Pauli operators are defined via X d = Ej=o L?® -0 1 
and Z d = YsjZo \j)( d Ul where Ctf is a rf -th root of unity 
rf = 1, and is addition modulo d. As for qubits, the 
qudit basis is = X d Z J d , with the commutation relation 
X° d Z k d = CdZdH and Xj = Zj = I. In the qudit case 
the Hadamard gate is replaced by a Fourier gate F, given 
by F\a) = d~ 1/2 E/c C fea \k), and satisfying F 4 = I and 
F 2 \a) = \— a) (—a is taken modulo d). 

Basic gates needed for our protocol are the qudit analogues 
of the CNOT (i.e. SUM), CPHASE and SWAP gates. Sup- 
pressing the dimension d and indicating the (control, target) 
qudits as subscripts |[T3l[T5l . 
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ilar in form to those found for the corresponding qubit gates: 
SWAP = CNOT 12 CNOT 21 CNOT 12 and CNOT 12 = 
(I ® H)CPHASE{\ ® H). 

Fig. [T] shows the simplified qudit SWAP gate when one of 
the inputs is |0). It builds an A^— qudit perfect mirror circuit 
A)] by replicating the simplified two qudit SWAP gate 
ljD) N — 1 times and adding redundant CPHASE and 
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Figure 1: Simplification of the qudit SWAP gate when one of the 
inputs is the |0) state. (A) Full SWAP from (kill, (B) shaded gates 
are redundant for input |0)„, (C) inserting FF^ 1 on either side of 
the target on CSUM gates, (D) simplification using ftc) . 



F gates, where 



F +1 |0). Through the repeated applica- 



tions of SF, with the global qudit gates F = Ylj=i F(j)< an< ^ 
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, one can transport the unknown quantum 



state from one end of the qudit chain to the other [21 1. As it 
stands, it would seem from Fig. |2jA) that perfect mirror inver- 
sion requires the particular initial state l^)®^)®^) • • -(g>|0). 
We now show that this is not the case. 

It was shown for the qubit case that the transport does 
not depend on the input states 0. The proof relies on 
showing that (HS) N+1 X l {a) = X l (N+1 _ a) (HS) N+1 and 

(HS) N+1 Z l {a) = Z l (N+l _ a) (HS) N+l which implies that 
any A^qubit density matrix is spatially inverted after the ac- 
tion of (H S) N+1 . For qudits, it suffices to show the analo- 
gous relations, i.e. 



F^"(F 1 S) N+1 X l f „^ = X, 
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We now develop the basic algebraic relations used to prove 
Q. These relations will be very important in a generalization 
from qudits to continuous variables, as showing that they also 
hold for the CV case immediately verifies the validity of the 
mirror circuit in that case. To simplify the notation, the over- 
bar denoting the global pulse is dropped unless stated other- 



wise by a subscript i.e. F 
calculation shows 



F, but F (a) -> F (a) . Direct 
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(F~ 1 S)X l (a) = X l {a _ 1) Z^ ) X l {a+1) (F~ 1 S) 

{F-'S)Z\ a) = X l ,s (F~ l S) . 
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To derive Eqs. |2} a graphical method |4| is used [Fig. 
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Figure 2: (A) Perfect mirror inversion circuit for qudits. Although 
particular input states are shown, the mirror action is independent of 
the input. (B) Evolution of the Z l (and X ) under repeated actions of 
F~ 1 S. The horizontal axis gives the number applications of F~ 1 S, 
while the vertical axis is the qubit position. The Z l (□) and X 1 (O) 
operators evolve according to (black shape indicating a negative 
power of the correponding operator). After N + 1 repetitions the Z l 
operator has moved to its mirror spatial position along the chain but 
has suffered I — > —I, Negative exponent is corrected by the final 



F operation. 



2[B)] to obtain (F^S^X^ = X^ +1 _ a) (F^S) 

\N+1 



\N+1 



and {F^S) N ^Z\ a) = Z^ +l _ a) {F^S)^\ It follows 

[via jyll] that extra final F ±2 corrects the exponent I — * —I. 
Fig. S[B) shows the evolution of an initial Z l operator after 
the consecutive action of F~ 1 S. It is extended to N + 2 so 
the evolution of both X 1 and Z l can be fully appreciated. It is 
important to note that the property Z d = I = X d was never 
used, thus allowing us to use the same method for the CV case 
which lacks this cyclic property. 

Continuous variable mirror transport:- For continuous 
variables ifTUl the same argument is used. The CV analogues 
of the generalised Pauli operators are X(q) = C,~ qp , Z(p) = 
( px where £ = e l / h , are non-commutative, X(q)Z(p) — 
q-ip Z(p)X(q), and have the following action on the compu- 
tational basis (position eigenstates x\q) — q \q))\ X(q) \s) = 
\s + q), Z(p) \s) — r] sp \s). By combining X(q) and Z(p), 
the displacement operator D(a) = exp(cra^ — a*a) can be 
formed. Given that any valid infinite dimensional density op- 
erator can be decomposed into a sum over coherent states via 
the P-function [16], p — J d 2 aP(a) \a) (a\, and that \a) = 
D(a) |0), the operators S(q,p) = X(q)Z(p); q,p g K 2 form 
a basis for Ti°° . This shows that CV- mirror transport protocol 
can be used to mirror any initial state of the CV-chain. 

Following the literature on CV-gates IfTUl . basic gates are 
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defined using the same notation as in qudit case. The posi- 
tion representation will be used for all operators. The Fourier 
transformation is F{a) \x) = (cr^/n) 1 J dy £ 2x y/< j2 sat- 



isfying F 4 = I and F 2 
basic gates needed are 



upon setting a = y/2. The 



u sum(12) = / dx \ x ) ( x \(i) ®X ( 2)(±x) = CN {12) (5a) 
Ucphase = dx \x) (x\ (1) (g) Z (2) {±x) = Sf 12} (5b) 



holds 



where the identity U CPHASE = F {2) CN {l2) F {2) 
again. The mirror circuit can then be built by using the + 
version of the defined gates and taking a = \/2. 

Note that using the + version of these gates allows to define 
the circuit independent of the basis [17]. With this in mind the 
reader shouldn't be surprised that the relations 

F ±2 Z(s)F ±2 = Z(-s), FZ{s) = X(-s)F, FX (a) = Z{s)F 

(6) 



F- l X{s) = Z(-s)F-\ F~ 1 Z(s) = X(s)F~ 



are also valid in this CV case. Since the algebraic rela- 
tions ( |4d| ), |4e} are the same, the same proof used for the qudit 
case is valid in the CV case. Alternatively, we can use the 
traditional representations of SUM and Fourier transform in 
CV IfTTl . i.e. without using the representation ( f5b]> ^ where 
the SUM gate can be represented as Usum(12) — CW(i2) = 
exp (— iqi ® p 2 ), which has the following actions: 

CN {12) : X x {x) ®l 2 ^ X l {x) ® X 2 (x) , (7a) 

Z x (p) -» Zi(p), X 2 (x)^X 2 (x) , (7b) 
U(g>Z 3 (p)^ Z- 1 (p)(g,Z,(p) , (7c) 
F : X —> Z, Z -> X- 1 . (Id) 

With these and F = Yif—i F(j), and S — Tl^i ^(jj+i)' we 
arrive at relations analogous to pdh and Bel F 1 S : X, 
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as desired. Thus ([8]l proves that the circuit shown in Fig. ^A), 
with the analogous CV gates performs perfect quantum mir- 
roring on a chain of CV systems. 

Circuit-QED quantum mirroring:- The above CV quantum 
mirror transport can be realised in Circuit-QED [12J. The goal 
is to spatially mirror quantum state of N—CV modes held 
in chain of superconducting coplanar waveguides (CPWs), 
nearest-neighbour-coupled by switchable Cooper Pair Boxes 
(CPBs). To this end, F and S ought to be executed, where for 
S, transformation cxp (— iq\ ® q 2 ) needs to be generated. We 
consider the arrangement where the CPB coupling adjacent 
1/4-wave CPWs can be switched between being biased by the 
difference or the sum of the voltages of the two CV modes. 
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Figure 3: Quantum evolution of the two CPW modes in the adiabatic 
limit for H-, taking \ip(r = 0)) = \a = 1, = 0, 1), (the two 
modes in coherent states \a) ® \/3) and CPB state is |1}), and taking 
w = 15GHz, Lu a/b = 3GHz, g a/b = 200MHz, 7 = 15MHz, 
K a/b — 1MHz. Plots show the reduced mode dynamics of the full 
master equation {9b} [Full], the effective master equation < | 1 0a| > [Eff], 
and subject only to the effective Hamiltonian {TOb} [Ham], with r = 
uiat/2-n, the nominal oscillation period of the CPW modes. 

The Hamiltonian and the resulting master equation for the 
two CV modes (operators a, b) and CPB (a± operators) oper- 
ating at the charge degeneracy point are ifTZll 



H±/h — u> a (ra + ujbb^b 



U 



g a (a + tf)±g b (b+tf) + (9a) 



P=j[H ± ,P]+ (£C(S-_) + fC(a) + fC(b)) p (9b) 

where H± is the dynamics for the case of sum or difference 
bias conditions, io a j h are the mode frequencies of the res- 
onators, g a j h are the couplings between the CPB and each 
CPW, luq is the CPB splitting, 7, K a /t are the decay rates for 
the CPB and two CPWs, and C(A)p = 2Apti - {A^A, p}. 

Assuming that CPB dynamics dominate all other time 
scales in (|9jl, i.e. ujq > LL> a /b,ga/b, an d 7 > K a /b, the fast 
dynamics of the CPB can be eliminated ifTSl . The equation 
for p = Ticpb{p}> describing dynamics of the two coupled 
CV modes is (for g a = 95) 

p = -i/h[H eff ,p}+ K j£(j)p + v£(s)p , (10a) 

j=a,b 

H e f f = Huj a a* a + fuo b Vb + hxs 2 , (10b) 

where s = X a ± X b and X a / b are the position operators 
of the two CV modes x — 5 2w o/[(7/2) 2 + ujq], and rj = 
(7 /2)/[(7 /2) 2 + ljq]. Fig. [3]shows that this approximation is 
valid for several hundred oscillator periods. 

To generate the unitary operator Uphase(12)> consider 
the effective Hamiltonian (lObi, with the mode frequencies 
and effective coupling strengths time dependent ui a (t) — 
oj b (t). Such rapid tuning of the resonator frequencies is 
achievable [19|. The canonical transformation H'(t) e ff = 
TH eff {t)f^ generated by f (a, b) = exp(-vr[at6 - Pa]/ 4) 
decouples ( |10b[ ) into two separate harmonic systems. The re- 
sulting Hamiltonian is H' eff (t) = u + (t)(a'^a' + b'^V) + 
2x(t)X 2 , i V , where u>+(t) = cu a (t) with either X 2 , for sum 
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exp(iX^, /10) for parameters 

99.9% 1 20 1 . Total time to effect desired unitary operator U' a is thus 
50 oscillator cycles. 



Figure 4: Pulse trains to effect U' a 
of Fig. ^to an accuracy of T — 



voltage biasing in (|9b| or X$, for difference voltage biasing. 



Under this canonical transformation the operator Uc phase 
also separates U' CPHASE = exp(iX^,) exp(— iX%,). Fur- 
ther, the Fourier gates applied homogeneously to both modes 
a and b transform to homogeneously applied Fourier gates 
on the modes a' and b . This gate is achieved via the natu- 
ral evolution of the uncoupled CPW resonators. We describe 
how to generate the U' CPHASE by showing how to generate 
the first component U' a = exp(iX^,): (i) numerically find 
time dependencies for w+(t) and x(*)> resulting in genera- 
tion of U' a while simultaneously b' acquires a known rota- 
tion exp(— iOb'^b); (ii) Switch to difference voltage bias- 
ing, evolving via a similar set of pulses, to effect U' b , while 
a! acquires an identical known rotation. These extra rotations 
are then absorbed into the preceding and succeeding Fourier 
gates in the overall application of the Fig.[2]A). To derive the 
required time dependent pulse we focus on the a! dynamics 
in H' e jj and move to a new time variable t = Qt, where 
w+(t) = w(l+Ci(t)/50),x(t) - x(l+C 2 (t)/50), and write 
H' a (r) - [1 + C 1 (r)/50]a'ta' + [ e + C 2 (r)/50]! Q 2 ,. C\{r) 
and Ca(r) are control functions and e = 2x/oj ~ 1 x 10 3 for 
the parameters of Fig. [3] To obtain pulse sequences which at- 
tain 99.9% accuracy in simulating U' a , we use GRAPE-based 
numerical optimization algorithms [20] in 50 cycles of the 
oscillator and with a Fock state truncation of 70 l22l . The 
pulses, displayed in Fig. [4] require tuning the CPW resonators 
by 2.4MHz and the CPB by 150MHz. It is clear that if one can 
generate larger detunings of the oscillators and CPB, the time 
required to accurately simulate U' a can be even shorter, thus 
allowing several mirror iterates to be executed before the res- 
onator decay degrades the dynamics. 

In summary a protocol for perfect quantum mirroring of a 
chain of coupled qudits or CV systems was proposed. Such 



capability will allow the quantum transport of entire entan- 
gled qudit or CV quantum states via quantum mirroring. The 
protocol uses only global applications of higher dimensional 
versions of the Hadamard and CPHASE gates. We anticipate 
that CV universal quantum computation can be achieved if 
the protocol is augumented with the capability of performing 
a conditional phase gate between the end two CV qubits in the 
chain. This would require a large cross-Kerr nonlinear inter- 
action between these CV systems. Finally, it was shown how 
to implement CV quantum mirroring in Circuit-QED. 
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